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Abstract 

When a foreign atom is placed on a surface of a metal, the surrounding sea of electrons responds 
screening the additional charge leading to oscillations or ripples. On surfaces, those electrons are 
sometimes confined to two-dimensional surface states, whose spin-degeneracy is lifted due to the 
Rashba effect arising from the spin-orbit interaction of electrons and the inversion asymmetric 
environment. It is believed that at least for a single adatom scanning tunneling microscopy mea- 
surements are insensitive to the Rashba splitting i.e. no signatures in the charge oscillations will 
be observed. Resting on scattering theory, we demonstrate that, if magnetic, one single adatom is 
enough to visualize the presence of the Rashba effect in terms of an induced spin-magnetization 
of the surrounding electrons exhibiting a twisted spin texture described as superposition of two 
skyrmionic waves of opposite chirality. 
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The lack of spatial inversion symmetry is the triggering ingredient for a wide range of 
new phenomena that are accessible with state of the art experimental techniques [THE]. 
Angle-resolved photoemission spectroscopy (ARPES) was the first tool used by LaShell 
and coworkers [7J to discover a small energy splitting in the sp surface state band of the 
Au(lll) surface. This splitting has been interpreted by the same authors as a realization 
of the interaction between the spin and orbital angular momentum, which promoted the 
renaissance of the Rashba physics. The surge of interest in similar effects involving the spin- 
orbit interaction in a structural asymmetric environment is the incitement of many additional 
sophisticated measurements and theoretical simulations on relevant surfaces with [HHID] and 
without [TTHT5] topologically protected surface states. Since charge oscillations (Friedel 
oscillations) induced by scattering of the surface-state electrons at a single adatom are blind 
with respect to such spin-orbit effects [16], an alternative [17] has been proposed on the basis 
of a multiple scattering study consisting on probing with the scanning tunneling microscope 
(STM) charge oscillations confined within a corral of adatoms. Theoretical calculations [IB] 
outline the possibility of very complex magnetic structures. We pursue a different route and 
investigate the possibility of grasping information on the spin-orbit interaction at surfaces 
exploiting the break of time inversion symmetry introduced by a magnetic adatom. As shown 
and discussed later, we discover an intriguing spin-texture in the spin-polarized electron gas 
surrounding the magnetic adatom (Fig. [TJ. This new magnetic behavior can be very similar 
to the topological twists, called skyrmions [IH3J I2II1 121] ■ Parts of our predictions have 
recently been verified in the interferences produced by the scattering at a MnPc molecule of 
the complex surface states of Bi(110) surface [T5] . 

Our investigation is based on a Rashba Hamiltonian [19] that describes a two-dimensional 
gas of free electrons confined in the (xy) plane of a metallic surface: 



considered with respect to a zero-energy reference defined by the bottom of the dispersion 
curve in absence of the spin-orbit coupling, m* is the effective mass of the electron and 
an is the effective Rashba parameter, describing the strength of the effect, whose value is 
determined in principle by the atomic spin-orbit strength as well as by the degree of asym- 
metry of the wavefunction imposed by the presence of the surface [16]. Here, however, an 
is chosen to model the experimental dispersion relation of the suface state. Quite gener- 
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FIG. 1: Skyrmionic-like spin-texture at the Fermi energy of Au(lll) surface, (a) Spatial visu- 
alization of the induced local density of states (LDOS) and the local magnetization direction of 
Au surface electrons after scattering with an Fe adatom. The spin texture found in (a) can be 
decomposed into a linear combination of two smoothly rotating skyrmionic magnetic waves with 
opposite vector chirality shown in (b). The wavelengths of LDOS oscillation, the left and right 
skyrmionic waves are ~ 18.7 A, ~ 17.3 A and 20.3 A, respectively. 

ally the gradient of the potential at the surface acts as an electric field E normal to the 
surface in the lab frame of the sample. Electrons propagating with momentum k across 
the surface experience this field in their local frame of reference as an effective magnetic 
field, B = -\-k x E, which is the origin of the functional form for the Rashba Hamiltonian 

m*c ' ° 

and defines a spin-quantization axis n(k) to be located in the surface plane normal to the 
wavevector k = (k x , k y ) = /c(cos</>, sin0) _L h(k), and c is the speed of light. The eigenvec- 
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tors | "01 (2)) °f the Hamiltonian (jlj), associated with the wavevectors k\ and fc 2 , are spin-up 
and spin-down states, respectively, with respect to n, but are also a coherent superposition 
of spin-up and -down states, (0i(2)(A;)) = e ^ (| t) — (+)«e i * > | 4-)), expressed by the spin func- 
tions | f) an d | 4), when measured with respect to the surface normal (z-direction). The 
eigenvalue spectrum E\(p)(k) = 2^(^1(2) — ^so) i s a two spin-split cone-shaped parabolic 
energy dispersion curve by which the origins of the parabola £1(2) are shifted with respect 
to k = by k so = ±m*afj/h 2 , i.e. k^ 2 ) = k + (—)k so (see Fig. [2|. When k 2 changes sign 
from positive to negative, the two branches of the dispersion curves cross. Unless stated 
otherwise, we shall consider in the upcoming text only positive k 2 . 

The scattering of the surface-state Rashba electrons on a magnetic adatom deposited on 
a substrate is investigated using scattering theory, which has been successfully applied in the 
description of electron scattering at adatoms [221424] . This theory involves the calculation 
of Green functions, which allow an elegant treatment of the electronic properties including 
the embedded adatom by solving the Dyson equation, G = G + G tG Q , where all quantities 
are site, spin and energy dependent matrices. We note that the electronic properties of 
the adatom are inscribed into the scattering matrix t, the amplitude of the electron-wave 
scattering at the adatom. Go is the Green function of the pure two-dimensional electron gas 
corresponding to the Hamiltonian (fTl) that is constructed using the eigenvectors \ipi(2))' 



Gq 



Gb e ^Gnd 



(2) 



In case the Rashba effect vanishes, the off-diagonal part of the Green function, Gnd, is 
zero and the diagonal part, Gd, reverts to the Green function of the free electron surface 
states. From G, quantities related to those measured by STM, such as the local density 
of states (LDOS) [25] n(r;E) = — 1/tt ImTr CT G(r, r; E) or measured by the spin-polarized 
STM (SP-STM) as the local magnetization density of states (LMDOS) [26] M(f; E) = 
— 1/n lmTT a aG(f, r; E) can be calculated, where a is the vector of Pauli matrices and Tr CT 
is a trace in spin space. It is convenient to express all quantities in cylindrical coordinates 
as it turns out that they depend only on the radial distance R measured from the position 
of the adatom. The LMDOS is expressed by two components, M(R) = (M z , M r ), the radial 
component, M r , and the one normal to the surface, M z . The azimuthal component, M^, 
vanishes. 
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FIG. 2: Energy dispersion of spin-split surface-state electrons. When plotted with respect to the 
two components of the two-dimensional wave vector k, the energy dispersion are cone-shaped. Due 
to spin-orbit interaction, two spin-split parabolas are obtained and are centered around ±fe so (see 
text). Arrows indicate the vector fields of the spin-quantization axes (or the patterns of the spin) 
at the constant energy contour. For every energy two opposite spins have different wave vectors 
leading to two concentric circle with opposite spins. The effective B-field felt by the electrons is 
perpendicular to the propagation direction defined by k. 

The obtained results are general, but to strengthen our point we consider as an application 
the system of an Fe adatom on the Au(lll) surface, since the system has been proven 
experimentally accessible by low-temperature STM and for Au the presence of a large Rashba 
effect has been shown experimentally [7J. On such a surface, Jamneala and coworkers [27J 
found that contrary to Ti, Co, and Ni, no Kondo peak was observed for V, Cr, Mn, and 
Fe adatoms meaning, that for the latter elements indeed a local magnetic moment exists 
and it is not screened by the conduction electrons [28]. We treat Fe as an adatom with 
a magnetic moment pointing along the z-direction perpendicular to the surface since our 
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first-principles calculations |29j (see supplementary material) predict this to be the easy axis 
with a magnetic anisotropy energy of about 12 meV. 

Hence, such a case imposes a scattering matrix that is diagonal in spin-space, 



(3) 

_o t u _ 

which can be related to the phase-shift 5(E) experienced by the incoming electron waves 
scattering at the adatom (t = h 2 /m*[exp(2i5(E)) — 1]). Note that a generalization to an 
arbitrary rotation angle of the magnetic moment will be straightforward using standard 
unitary transformations. For Fe, as for Au, all majority-spin states are fully occupied and 
the scattering of majority electrons in the vicinity of the Fermi energy (Ep) is practically 
zero, i.e. the phase-shift vanishes (t^ = 0). The minority-spin LDOS shows, on the contrary 
to Au, high values around the Fermi energy and therefore a large scattering. Thus, we 
assume a phase-shift of n/2 (t^ = —2h 2 /m*). 

The parameters describing the Au(lll) surface state are identical to those chosen by 
Walls and Heller [T7j: m* = 0.26 m e , Ep = 0.41 eV and an = 0.4 eVA that correspond to a 
Fermi wavelength Ap = % ~ 37.4 A and a spin rotation length A so = -r- of 230.5 A. 

After solving the Dyson equation, the induced circular LDOS oscillations, An(R; E), em- 
anating from the adatom located at the center (R = 0) is given by An = — 1/n Iul[(G'dG'd + 
G r NDG r ND)(^tt + ^44-)] an d is plotted in Fig. [IJa) for all states at E = Ep. At large dis- 
tances R, the circular induced standing wave undulations, An(R;E), can be expressed as 
~ lz 2^ k 2 R \/k\k2 cos(2kR). When the spin-orbit coupling is negligible, i.e. k\ ~ k 2 , we recover 
the conventional form of the adatom induced energy dependent charge density oscillations 
~ c °2*fcj^ [221 123] . Regardless of the spin-orbit interaction, there is only one wave vector, 
2k = k\ + &2, that describes the oscillations, with the corresponding wavelength at Ep given 
by \pj2 ~ 18.7 A . Consistent to the work of Petersen and Hedegard [16] and confirmed 
by Walls and Heller [IT] , no signature or information on the Rashba effect or spin-orbit 
interaction of the surface atoms can be extracted from the LDOS. 

This picture changes fundamentally when we look at the LMDOS. Due to the presence of 
a magnetic adatom [t^ ^ t^), the magnetization density perpendicular to the surface and 
in the surface plane are non-zero. We find concentric rings M(R) of equal size magnetization 
density surrounding the Fe adatom, with magnetization densities emanated at the center of 
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the Fe atom wrane ling in the (M Z (R), M r (R)) plane (Fig. [LJa)). Contrary to the LDOS, 
the magnetization density is non-trivially modulated by the Rashba effect. Analysing the 
asymptotic behavior of M Z (R) 

~ 2 7r 2^2fc2^ (^i cos(2,k\R^j -\- k<i cos(2/u2-/?) ) as calculated from 
M z = -l/nlm[(G D G D -G ND G ND )(t n -t u )]) as well as of M r {R) ~ ^^{h sm{2hR)- 
/c2 sin(2/c2-R)) as calculated from M r (R) = — 2/n Im[GD(/tt ~~ ^u^nd], one finds that both 
wave vectors k\ and /C2 enter nontrivially. This observation implies the possibility of using 
a SP-STM with an out-of-plane magnetized tip to probe the ^-projected Rashba induced 
interferences. Any magnetic signal detected by a SP-STM with an in-plane magnetized tip 
would be a clear fingerprint of the Rashba effect. If there is no Rashba effect, i.e. G*nd = 0, 
the quantization axis is only determined by the adatom and the scattering electrons are not 
a coherent superposition of spin-up and down-electrons that that led to a finite in-plane 
spin-component. 

Such spin-textures are a reminder of skyrmionic magnetic configurations but are com- 
pletely different. Indeed, with skyrmions, the magnetic moments can rotate smoothly from 
one direction at the center of the structure to the opposite direction. Here, however, the 
magnetization ripples experience phase-switch as well as beating as sketched in the top- 
right of Fig. [3j In order to simplify our analysis, the induced LDOS and the z- and 
radial components of the magnetization passing by the adatom position are additionally 
illustrated in the top-left of the same figure. Surprisingly, the amplitude of the M r and 
M z oscillations are of the same size with a phase-shift between them that depends on 
the strength of the Rashba coupling term. This can be better understood when looking 
at the asymptotic behavior. For instance, M z can be rewritten in terms of k so and k as 
~ n 2fc> k 2 R (k cos(2kR) cos(2k so R) — k so sm(2kR) sin(2k so R)). At the Fermi energy of gold, 
k so is expected to be very small compared to k. At distances much smaller than the spin- 
rotation length but long enough that the previous asymptotic behavior holds, M z simplifies 
more to w 2™2 kR cos(2kR) and oscillates similarly to the LDOS with or without spin-orbit 
interaction. This explains the wavelengths of the induced initial ripples being, as expected, 
twice the Fermi wavevector. In fact one has to go far beyond the vicinity of the adatom, at 
distances of the order of the spin-rotation length, to see an effect achieved by the intriguing 
phase switch observed around 60 A. Indeed, the new spin-rotation length corresponding to 
a full rotation from 0° to 180° occurs around 115 A, which means that the change of sign 
occurs at half this value. Another substrate with a stronger Rashba coupling parameter 
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would impose a more important phase-shift at smaller distances. M r , which is proportional 
to — ^5ppfi (k cos(2fc-R) sin(2fc so -R) + k so sm2kRcos(2k so R)) can be further simplified - for 
rather short distances - to — ^tj^f - cos(2kR), similar to the cosine behavior that characterizes 
M z . 

Although the spin-texture looks very complex at first sight it can be understood 
as a linear combination of two skyrmionic waves, M = M kl + Mk 2 , with M^^R) oc 
ki(cos(2kiR) , — sia(2kiR)) and M k2 (R) oc fc 2 (cos(2/c2-R), sin(2fc 2 -R)) of opposite vector chi- 
rality defined as c = M(R) x M(R + dR). The chirality c k2 = sin(2k 2 dR)e ( j ) points in the 
0-direction and forms with the directions e z and e r a right handed coordinate system, while 
Cfcj = — sin(2kidR)e r f ) points in the — 0-direction and forms with e z and e r a left-handed 
coordinate system. Accordingly, we call M kl and M k2 a left- and right-rotating skyrmionic 
wave, respectively. Our definition of a skyrmion in the actual work designates well-defined 
magnetic waves which are (i) centered around the adatoms, and have (ii) fixed rotation 
sense [Ij. One great virtue of skyrmions is their topological nature that is protected under 
the presence of reasonable size magnetic fields. 

This multi-skyrmionic magnetic waves can be destroyed or manipulated by modifying 
the spin-orbit interaction. By switching it off, only the magnetization along the ^-direction 
would be finite. One way to tune the Rashba effect is to change the substrate nature 
gradually [301 EE]- Such an experiment is very difficult to realize but we propose another 
tuning method: Instead of the Fermi energy, one could probe at different energy varying the 
bias voltage between tip and sample. For example, by decreasing the value of the energy 
probed, the wavelength corresponding to 2k decreases which is accompanied by a decrease 
in the number of nodes. This is observable in the two additional examples depicted in 
Fig. [3] calculated at 140 meV and 20 meV. The final multi-skyrmionic waves have a different 
texture but share the common feature in the phase switch at 60 A. 

One should bear in mind that, as for the LDOS, two different regimes have to be expected 
for the magnetization components depending on the sign of k% (see supplementary material). 
This suggests that the scanning tunneling spectroscopy (STS) experiment of Ast et al. [13] . 
but spin-polarized, would be also a way to verify our predictions. 

To summarize, we have revealed a new type of spin-texture (see Fig. [4]) induced in a 
confined electron gas subject to the Rashba effect by a magnetic atom. These structures 
can be understood as a kind of combination of skyrmionic-like waves of opposite chirality. 



8 




(I 211 40 till lUl 



Distance from the adatom (A) 



FIG. 3: Comparison of the induced LDOS and LMDOS components. Due to the Rashba effect 
the Fe magnetic adatom induces non-trivial spin interferences: on the left panels are shown the 
radial dependence of the LDOS and LMDOS at different energies while on the right panels the 
corresponding magnetization unit vectors are depicted. 

By tuning the energy level or the spin-orbit strength, the observation of the final spin- 
texture is possible because of the large magnitude of the induced magnetization and the 
creation and exploration of new complex configurations that manageable with state of the art 
experiments. Such non-trivial magnetic Friedel oscillations have an impact on the magnetic 
interactions between adatoms and other nanostructures and consequently on their magnetic 
behavior. Simple asymptotic expressions for the induced magnetizations are derived which 
offer a simple understanding of new types of experiments involving the manipulation and 
construction of adatom based spin-nanostructures. 
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and Drs. P. Lazic and H. Schumacher for their assistance in preparing some of the figures. SL 
acknowledges the support of the HGF-YIG Programme VH-NG-717 (Functional nanoscale 
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